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A single superconducting artificial atom can be used for coupling electromagnetic fields up to the single pho¬ 
ton level due to easily achieved strong coupling regime. Bringing a pair of harmonic oscillators into resonance 
with the transitions of a three-level atom converts atomic spontaneous processes into correlated emission dy¬ 
namics. We present experimental demonstration of two-mode correlated emission lasing in harmonic oscillators 
coupled via a fully controllable three-level superconducting quantum system (artificial atom). The correlation 
of emissions with two different colors reveals itself as equally narrowed linewiths and quenching of their mutual 
phase diffusion. The mutual linewidth is more than four orders of magnitude narrower than the Schawlow- 
Townes limit. The interference between the different color lasing fields demonstrates that the two-mode fields 
are strongly correlated. 

PACS numbers: 42.50.Lc, 42.65.Lm, 03.67.-a, 85.25.Cp 


Excited free atoms relax to their ground states via incoher¬ 
ent spontaneous emission processes [1]. In an ensemble of 
three-level atoms - the fundamental objects of quantum op¬ 
tics - coupling of the lowest pair of levels to a harmonic os¬ 
cillator leads to conventional lasing if incoherent relaxation 
from the second excited state to the first excited state creates 
population inversion. The situation is dramatically changed 
when both transitions of the atoms are coupled to indepen¬ 
dent oscillators: the spontaneous emission processes are re¬ 
placed by coherent energy exchange and correlated dynamics 
between the oscillators, known as correlated emission lasing 
(CEL) [1-6]. The correlated dynamics reveals itself as a re¬ 
duction in relative random phase-diffusion noise, resulting in 
the suppression of mutual peak width below the Schawlow- 
Townes limit [7]. The suppression of phase-diffusion noise 
has been experimentally observed for two polarization modes 
in a HeNe laser [4-6]. However, the best result observed so 
far has been about 2.3% of the Schawlow-Townes limit for 
an ensemble of natural atoms [6]. Recently, the classical las¬ 
ing effect has been reproduced on single atoms [8] or single 
superconducting quantum systems (artificial atoms) [9-11], 
in which strong coupling to circuit elements can be easily 
achieved [9-13]. In this Letter, we demonstrate the coherent 
dynamics of two harmonic modes with different frequencies 
in a transmission line resonator (TLR) coupled through a sin¬ 
gle three-level artificial atom. Oppositely to the ensemble of 
atoms with uncorrelated dephasing in different ones, the sin¬ 
gle atom allows to ideally suppress the phase diffusion. We 
demonstrate quenching of the mutual phase diffusion to the 
level better than 10 -4 . The correlation between two lasing 
fields is displayed in the interference between them. 

In addition, we point out that our system presents a fun¬ 
damentally new circuit. In usual quantum electrodynamics, 
atoms are coupled through a resonator mode, however the re¬ 


versed circuit with linear oscillators coupled via transitions of 
a quantum system is very difficult, if not impossible, to real¬ 
ize with natural atoms, spins or quantum dots. Such a system 
presents a novel approach and allows to demonstrate a series 
of qualitatively new phenomena of quantum optics. 

As shown in Fig. 1(a), Fig. 1(b) and Fig. 1(c), the single 
three-level artificial atom with cyclic transitions [14] is based 
on a “tunable gap flux qubit” circuit [15, 16] capacitively cou¬ 
pled to a multimode transmission line resonator. Such a cou¬ 
pling for the flux qubits has been realised in [17] and the tran¬ 
sition matrix elements for the three-level system calculated in 
[18] are consistent with our experiment. The artificial atom is 
situated in the voltage antinode of the resonator. The transi¬ 
tion frequencies (c o eg , ujdg and ujde) between the three lowest 
levels of the atom | g), \e) and | d) (the ground, first excited 
and second excited states) are controlled by an external mag¬ 
netic flux, reaching minimal uj eg and uodg at half-integer flux 
quanta 4 >at = (N + (where N is an integer). Compared 
with systems based on the conventional flux qubit geometry 
[19, 20], our atom has additional tunability owing to the im¬ 
plementation of an a-loop - a dc SQUID, which allows the 
flux tunnelling energy A to be changed by choosing By 
selecting 4 >at and additionally tuning 5<f> (<C 4>o)> one can ad¬ 
just the transition frequencies uo eg and ujde to the resonance 
with the two lowest modes of the resonator - the two indepen¬ 
dent fixed-frequency oscillators r\ and r 2 with the frequencies 
(jj\ and cc 2 , respectively. 

The entire system presented in Fig. 1(d) can be described by 
the state \N\nJNf) = | Nf) ® \n) 0 | N 2 ) with the three quan¬ 
tum numbers n, Ni and 7V 2 , denoting the three-level atomic 
states | n) = {\g),\e),\d)} and the Fock states of resonators 
|7 Vi j2 ) = {|0>, |1), |2),...} with the photon occupation 7 Vi j2 
in the first and second oscillators. The full Hamiltonian of the 
two harmonic oscillators coupled via transitions of the three- 
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FIG. 1: Device description, (a) Optical micrograph of the device. 
The meandering structure is a TLR made of a 50-nm-thick Nb film 
capacitively coupled to external coplanar waveguides, (b) Magnified 
micrograph of an artificial atom based on the tunable gap supercon¬ 
ducting flux qubit geometry. The o-loop (with red junctions) allows 
the flux tunnelling energy to be tuned, (c) Simplified circuit model 
of the atom capacitively coupled to the TLR through capacitance C c . 
The atom is effectively coupled to two oscillators with the frequen¬ 
cies cji and u 2 , schematically shown by the dashed lines, (d) Prin¬ 
ciples of the device operation. The transitions |cZ)-|e) and |e)-| g) 
are resonant with the first and second oscillators, respectively. The 
transition \g)-\d) is driven by a classical microwave pumping field. 


level atom and driven by a microwave with the frequency 
uj p « c j dg and amplitude 2 HQ is 


H = huia\ai ± /Lu 2 <4 a 2 ± kw dg a dd ± huj eg a ee (1) 
±%1 (a\&ed ± Cde'j ± %2 a ge + a 2& eg^ 

+ 2 Ml((J d g + (Jgd) COS UJpt. 

Here, <jjk = \j)(k\ is the projection/transition operator for 
the atomic states with {j, k} = {e, g , d}. a\ (ai) and a\ (a 2 ) 
are the creation (annihilation) operators of the first and second 
oscillators, respectively. The atomic states | g) and | d) are cou¬ 
pled by the pumping field. To eliminate the time dependence, 
the Hamiltonian can be further simplified using the rotating 
wave approximation, in which the two oscillators appear to be 
resonantly coupled via the quantum system. The dissipative 
dynamics of this artificial atom-resonator coupled system is 
described by the Markovian master equation with accounting 
energy relaxation and decoherence in the system [18]. 

We perform microwave characterization (transmission and 


emission measurements) of the device in a dilution refriger¬ 
ator at a temperature of 30 mK. The device fabrication pro¬ 
cess and experimental setup are described in detail in Sup¬ 
plementary Materials [18]. From the transmission measure¬ 
ment, cci/27t « 6.0016 GHz and cj 2 /27t « 11.9979 GHz with 
the corresponding linewidths (decay rates) in the oscillators 
of tti/27r = 0.63 MHz and ft 2 /27r = 1.94 MHz at the base 
temperature. We find that 4>jv ~ 1.54> 0 ( A = h x 1.51 GHz) 
and 5$ ~ 18 x 10 -3 <f>o is one of the best choices to reach 
the desired double resonance. To characterise the atom- 
resonator interaction, we measure the normalized transmis¬ 
sion through the system \t/to\ around uj\ and cc 2 as a func¬ 
tion of the probing frequency and flux bias S 4> using a vec¬ 
tor network analyser (Fig. 2(a) and Fig. 2(b)), where to is the 
maximal transmission amplitude. Anticrossings due to vac¬ 
uum Rabi splittings are clearly observed when u eg « in 
Fig. 2(a) and u eg « cc 2 in Fig. 2(b), respectively. The cou¬ 
pling strengths of go = 2tt x 36 MHz (corresponding to the in¬ 
teraction |0e0) <=> \lg0) at « ±9 x 10 -3 <f> 0 ) and g 2 = 
2t r x 78MHz (|0e0) <=> \0gl) at S <f> « ±19.5 x Hr 3 4> 0 ) 
are obtained by fitting the anticrossings. 

In Fig. 2(c), we use the two-tone microwave spectroscopy 
technique to find the energy levels of the atom-resonator cou¬ 
pled system [21]. We measure the transmission at a fixed 
frequency close to and sweep the pumping frequency in 
the range from 1 to 16 GHz (limited by the bandwidth of the 
entire system). However, we show our theoretical calcula¬ 
tions (which are in a good agreement with the experimen¬ 
tally measured lines) above the bandwidth limit by dashed 
lines. The applied probing power is sufficiently low to keep 
the average number of photons inside the resonator less than 
one ((TVi) < 1). At 5§ / 0, all the transitions between 
the three levels are allowed [14]. The black dashed lines in 
Fig. 2(a), Fig. 2(b) and Fig. 2(c) are the fitting curves obtained 
from the Hamiltonian of the atom-resonator system [17]. The 
dashed red box in Fig. 2(c) shows the avoided crossing due 
to the interaction between \0d0) and |0el) characterised by 
#4 = 2tt x 205 MHz at ~ ±2.5 x 10 _3 4> 0 , where the 
second excited atomic state is converted into a photon in r 2 
and the first excited atomic state. Note that g 4 extracted from 
the experiment is also consistent with the theoretical calcula¬ 
tions [18]. 

To couple the oscillators via the atomic transitions, we set 
the bias to ^^>5 = — 18 x 10 _3< F 0 (red dashed line), where 
uj eg ^ uo 2 and ujde ~ cji, and apply the external microwave 
drive at a pumping frequency of uo p = 2tt x 18.0055 GHz. 
Although gi cannot be directly measured, we take its value 
as 27r x 90 MHz from simulations, which well reproduces 
other couplings (detailed in Supplementary Materials [18]). 
Strong emission peaks appear when the pumping power ex¬ 
ceeds P « —100 dBm. The emission disappears, when the 
atom is detuned from the double resonance. We assume that 
the threshold has a different nature from conventional lasing 
and is due to the remaining detuning between cc^ e , uji and/or 
cjJ eg , cc 2 , compensated by Rabi splitting, which according to 
our estimates is 11/27r ^ 900 MHz at Ud g - This is supported 
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FIG. 2: Energy levels of the three-level atom-resonator coupled system, (a) and (b) Normalized transmission amplitude through the resonator 
modes as a function of 5&. The black dashed lines are theoretical fittings from the full system Hamiltonian, (c) Intensity plot of the phase of a 
transmission signal as a function of driving frequency and 5&. The black dashed lines represent theoretical calculations. Our goal is to study 
the double resonance shown in Fig. 1(d), which takes place under the conditions marked by the red squares. 


by our simulations [18]. Note also that we expect another 
resonance close (but not equal) to uj p . However it does not 
affect the overall picture because the drive is essentially clas¬ 
sical with many photons in the mode. The emission spectra 
of the coupled system in the vicinity of cci and CC 2 are simul¬ 
taneously monitored by two spectral analysers (SAs). We use 
the techniques developed in [22] to eliminate the noise back¬ 
ground of the low temperature amplifier by subtracting traces 
of OFF from ON emission spectra. The results are shown in 
Fig. 3(a) (black open triangles) and Fig. 3(b) (black open cir¬ 
cles). Importantly, the peaks are well fitted by equal width 
Lorentzian curves with Acc e i « Acc e 2 = 27r x 0.80 MHz, 
indicating the interaction between the oscillators. Note that 
Acc e i 5 2 is about two times lower than ft 2 and three times 
lower than the effective relaxation rate of the entire system 
of fti + ft 2 ~ 27r x 2.57 MHz. The center frequencies of the 
two emission spectra of uj e i = 2 ttx 6.00086 GHz and cc e 2 = 
27tx 12.00471 GHz are somewhat shifted from uji and UJ 2 by 
uj e \ —uj\ = —27r x 0.74 MHz andcc e 2 — CC 2 = 27r x 6.81 MHz, 
which is probably due to dispersive shifts from the detuning 
between the artificial atom and the resonator modes [17]. We 
emphasise that cc e i + ^e 2 = w p holds with high accuracy. 
Although the emission frequencies slightly depend on other 
parameters, such as pumping power and frequency, the above 
equation still holds [18]. 

The estimated total emission powers, obtained from the ex¬ 
tracted area of the emission curves, are Pi « —134.4 dBm 
and P 2 ~ —129.4 dBm, with an accuracy of about 3 dBm due 


to the uncertainty in the calibration of our setup. They roughly 
correspond to (Ni) « 5 photons and (7V 2 ) ~ 2 photons at r\ 
and 7 * 2 , respectively. (The photon number is derived using the 
equation (Nk) = 2P/ c /(fkj/ c ft/ c ), where k = {1,2} and the 
prefactor 2 is due to the escape of photons from either end of 
the resonator with equal probability.) The ratio (Ni)/(N 2 ) is 
found to be close to the expected value as the photon creation 
rate in each mode under stationary conditions is expected to 
be given by (Ni)ki = (7V 2 )ft 2 , which originates from the 
energy conservation law: each absorbed photon at uo p is split 
into a pair of photons at uj e 1 and cc e 2 . 

To demonstrate quenching of the phase diffusion in the two¬ 
mode CEL, we mix the amplified emission signals at cc e i and 
cc e 2 using a conventional microwave mixer. Then, we mea¬ 
sure the signal around the sum of the frequencies cc e i + ^e 2 
by the SA [18]. Note that the phase diffusion in our case is 
different from the relative phase diffusion in Refs. [4-6] and 
we call it mutual phase diffusion. To avoid direct leakage of 
the pumping tone through our circuit, we filter out high fre¬ 
quencies before mixing and confirm the absence of any out- 
coming signal at uj p . Next, we simultaneously monitor the 
signals at cc e 2 (cc e 1 ) and uj e 1 + w e 2 - The result shown in 
Fig. 3(c) demonstrates a strong narrow sum signal at a fre¬ 
quency of exactly 00 p only when emissions at cc e i ,2 occur. De¬ 
spite mixing a pair of signals of Acc e width, the resulting peak 
width of Acc sum = 2tt x 9.4 Hz is found to be limited by the 
bandwidth of the SA (10 Hz). We note that the Schawlow- 
Townes condition for the linewidth is expected to be modi- 
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FIG. 3: Emission spectra of two-mode CEL and quenching of the 
phase diffusion noise, (a) and (b) Normalized transmission through 
oscillators 1 and 2 (red closed and open rhombus) and emission 
spectra from the oscillators (black open triangulars and circles). 
The emission spectra have the same widths of Ac<7ei = Acj e 2 = 
27r x 0.80 MHz. (c) Quenching of the mutual phase diffusion. The 
two frequency emissions are amplified and mixed. The resulting up- 
converted signal (uncalibrated) is very narrow and limited by the SA 
bandwidth. 



FIG. 4: Correlation between the two emission modes, (a) Construc¬ 
tive (red dots) and destructive (black dots) interference peaks of the 
emissions with two different colors. To observe the interference, the 
uo e 2 -signal is downconverted to uj p — uo e 2 and a variable delay is in¬ 
troduced in the uj e 1 emission signal before summing it with u p —uo e 2 - 
(b) Interference fringes of the peak amplitudes from (a) as a function 
of the phase delay between the two signals (red dots). The ampli¬ 
tudes of both signals are adjusted to obtain nearly maximal modula¬ 
tion. The intensity derived from the experimental peaks (red dots) is 
fitted by a sine function (black curve). 


fied in our system to Accst = (ft 1 + ft 2 )/( 2 iVt 0 t), which is 
about 250 kHz for 10 photons. The experimentally measured 
linewidths of different color emissions satisfy the condition 
fti + ft 2 > Acc e i ,2 > Accst, which usually takes place for 
conventional lasers due to additional broadening mechanisms. 
However, the mutual phase diffusion noise (characterised by 
Acj sum ) is quenched being at least four orders of magnitude 
lower than Ao;st, which can be explained by the theoretical 
model [18]. 

To further quantify the correlation between the emissions, 
we downconvert the emission at cc e 2 by mixing it with the uo p 
tone and select out the difference frequency uo p — cc e 2 . Then, 
we linearly add the resulting signal to that at uj e 1 using a lin¬ 
ear adder and measure the resulting signal by the SA [18]. The 
intensity of the observed signal is extracted by fitting with a 


Lorentzian curve (red curve in Fig. 4(a)) and then is studied 
as a function of the phase delay introduced in the cj e i emitted 
signal. As shown in Fig. 4(b), we observe the sinusoidal os¬ 
cillations of the intensity of the interfering signal as a function 
of the delay, indicating a strong correlation between the two 
lasing fields. The interference is an additional evidence for 
the quenching of the mutual phase diffusion. 

Finally, we briefly summarise other aspects of the system, 
consisting of the two oscillators coupled through the transi¬ 
tions of the three-level atom, which can be verified experi¬ 
mentally in the future. The single three-level atom CEL sys¬ 
tem provides an ideal testbed for studying the nonlinear cou¬ 
pling of light in cavity modes. It was theoretically found that 
in contrast to two-mode squeezed states by parametric am¬ 
plifiers [23-27], a bright and arbitrarily high degree squeezed 
lights can be generated inside a cavity by a system of atoms in¬ 
teracting with a two-mode cavity because of the strong atomic 
nonlinearity [3, 28, 29]. Moreover, the nonlinear coupling of 
light in the two cavity modes can lead to two-mode squeezed 
states [1]. We calculate the time dependence of the sum of the 
quantum fluctuations of the two lasing fields (A u) 2 + (Av) 2 
based on the experimental parameters and find that it can be 
smaller than two [30-33], which suggests the possible exis¬ 
tence of quantum correlation between the two fields in our 
sample for the short time after tuning on the pump (details are 
given in Supplementary Materials [18]). In addition, the quan¬ 
tum coupler provides an ideal (the strongest possible) nonlin¬ 
earity for weak fields up to the single quantum (photon) level 
in the strong coupling regime. For example, there is a the¬ 
oretical discussion about strong coupling of a phonon mode 
and a photon mode with a single superconducting three-level 
artificial atom at single-photon level [34]. 

In conclusion, we have demonstrated correlated dynamics 
between two oscillators coupled via transitions of a three- 
level macroscopic artificial atom. The dynamics is observed 
as the quantum noise quenching of the mutual phase diffu¬ 
sion in the correlated emission lasing. Our results are the first 
step to generate the quantum correlation and entanglement be¬ 
tween two lasing fields with a single atom, which is an impor¬ 
tant step towards constructing a quantum network in the mi¬ 
crowave domain. This work may also lead to an improvement 
in ultrahigh-sensitive interferometric measurements^]. 
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TWO-MODE CEL THEORY 


the cavity field 


Master equation of reduced density matrix for two cavity modes 

In the rotating reference frame with U = exp(— iuj p addt — 
iujpa+ait), the Hamiltonian for the three-level system inter¬ 
acting with two cavity fields (Eq. (1)) is simplified to 


H = HAia[ai + h(jJ 2 a\a 2 + kdjJ eg G ee # hA 2 G d d 

+ h[flexp(i(f))(Tgd + gia\cr ed + g 2 alG ge \ + H.c(Sl) 


with two detunings, Ai = and A 2 = uJdg — w p . 

Taking the reduced density matrix of the three-level system 
and two cavity fields, p, we can write the master equation 


dp 

dt 


~j^[H)P\ H — [2(7gdP^dg t7dg®gdP P&dg&gd] 

732 

H [2(7 e dP(7de tTde^edP P&de<7ed\ 

^21 

H Zff [2(7g e p(7 eg (7 e g(7gep P&eg® ge\ 

i ^ 22 ro i 

+ ^ [^eeP&ee (7eeP P&eel 

733 

+ — [2(7ddp(7dd - (7ddP - p(7dd\ 

— a\aip — pa\ai^ 

+^ 2 ^2a 2 pa2 - a\a 2 p - pa\a^j , 


-\-Ki (^ 2 aipa\ 


dpj_ 

dt 


+ 

+ 

+ 

+ 

# 

+ 


—i[Aia\ai + uj 2 a\a 2 , py] (S2) 

{ a lP ( gg,A + r*lPdg,A ( a 2P/<4 “ Pf a 2 a 2) + H.C.} 

{ [“2 Pm,A + M 2 P { gd,A ( a \pf a i - Pfaia\) + H.c.} 

{ a iP ( eeN a \pf a 2 - a 2 a\p f ) + H.c.| 

( a 2/4°+ a iP/ a i - a\a lPf ) + H.c.j 

{ (fUf d ,A + v *2Pm,a) (« 2 p/oi - p/aia 2 ) + H.c.} 

{( 'PzP^Ia + V 1 ++) ( a 2 P/«i - axa 2 pf) + H.C.} 

{^ 2 PeJ A (alp f a\ - p f a[al) + H.c.} 

{viP { el A (. a iPf a 2 - Pfaia 2 ) + H.c.} 

2 

Y Ki ( 2a ^/ a ! - a l a iPf - Pf a \ a i)- 


i=1 


as 




(0) 

P 9 Ia 


(731 + 732)721(^2 + ^r 2 ) + r 2 7 2 i^ 2 
(731 +732)721(^2 + |r 2 ) + ^ 2 r 2 (27 2 i +732)’ 

_ 732 ^ 2 r 2 _ 

(731 +732)721(^2 + \ F 2 ) + ^ 2 r 2 (27 21 + 732)’ 
(731 + 732)721 ( iA 2 + lr 2 )U 2 exp(i 0 ) 

(731 + 732 ) 721(^2 + |r|) + ^ 2 r 2 (2 72 i + 7 32 ) ’ 


where the decay rates are 


with the decay rates k>i and k 2 of the two cavity modes and 
7 jk of the atom, where {k,j} = 1,2, 3. Using standard meth¬ 
ods of laser theory as in Ref. [1], we can eliminate the vari¬ 
ables of the three-level system and keep to the first order of the 
coupling constants gi (i = 1,2). The following is the master 
equation for the reduced density matrix of the two modes of 


Ti = 721 + 722 , 

1^2 = 733 + 732 + 731, 

r 3 = 733 + 722 + 721 +732 + 731- 
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O*- Also ol\, fa, i^i and Vi with i = 1, 2 are given by 

a i = ^A 3 + f (S3) 

<^2 = ^—^21 + — iAi^ , (S4) 

Pi = ^ ()a 3 ! ^ iA^ , (S5) 

P‘2 = ^-iW2 + 21 + y) , (S6) 

and 


Mi 


(S7) 

M2 
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= iSl-±exp(i<f>), 

U 2 

(S8) 


= D * ex p(^)> 

(S9) 

^2 

= D * ex p( *<A)- 

(S10) 


The denominators D\ and £>2 in the above equations are 


D 1 


D 2 


iw 2 - iu 21 + tt 


A ^ 
-2CJ21 - «Ai + — 


+ iA 3 H——^ + ^ 2 , 

^ (iA 3 — iAi + + ^ 2 , 


Eqs. (S11) and (S12), respectively. We then can also obtain 

= —2Ki(a\ai) + {^2 pfl A + «2 Pm,a + Hc } ( S13 > 
+ { (l4p ( gd,A + a 2 Pm,a ~ ^P^a) <4 a i> + Hc } 
+ { [ P lP ( gg,A - V lp { ee,A + P^P^a] («1«2) + Hc } > 

= — 2ft 2 (4 a 2) 3“ (^1 *H“ Qfi )Pee,A (S14) 

+ { [«1 P { ee - «1 P { gg ~ Pipfd,A (4 a 2) + H.C.j 

- { [AmSa + - ^pSa] ( a l a 2) + H - C -} * 


d(<4a 2 ) 

dt 


= i(cji + cu 2 -u d )(a\al) - (k x + « 2 )(44) ( S15 ) 
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The steady-state average photon number can be approxi¬ 
mately obtained from the above equations. At the working 
point, we expect A 2 = 2tt x 30 MHz, U = 2tt x 900 MHz, 
0 = 0, cJ eg = 2ttx 11.4979 GHz andcc^e = 2ttx 6.5376 GHz. 
Using the dissipation terms of the two modes and the artificial 
atom (the dissipation terms of the atom are extracted below), 
we obtain (a^ai) ~ 5 and (a 2 "a 2 ) ~ 2. Furthermore, both 
(afai) and (ajc^) decrease when we change the pumping 
power or pumping frequency away from these optimal param¬ 
eters, in good agreement with the experiment. 


with the detuning A 3 = Ud g — u p — w e g- Using the mas¬ 
ter equation for the reduced density matrix of the two cavity 
modes, we have 


d{ai) 

dt 


d(a 2 ) 

dt 


-iAi(ai) - Ki(ai) + j±-a 2 [p^U - p { °J iA ]( a i) 
P*2p ( gd,A( a l) - V l\p ( ee,A ~ P^A^ 

^2P£U<4>, (SH) 

—iui 2 {a 2 ) - K 2 (a 2 ) + ^al\p^A ~ CaK^) 

n rt (°) In \ 5l52o*.(0) , t\ 

PiPgd,A\ a 2 ) jj* PlPdg,A\ a l / 

^2 [p^A ~ Pm A (4) ■■ ( S12 ) 


Quantum correlation between two lasing modes 

We use a criterion derived in [2] to study the quantum cor¬ 
relation between the two lasing modes in two-mode correlated 
emission lasing. According to this criterion, a necessary and 
sufficient condition for a system to be in a quantum entangle¬ 
ment state is for the sum of the quantum fluctuations of the 
two Einstein-Podolsky-Rosen (EPR)-like operators u and v of 
the two modes to satisfy the inequality 

(A u) 2 + (At)) 2 < 2. (S16) 

Here, u = x\ + x 2 , v = Pi — P 2 , %j = (% + a+)/\/2 and 
Pj = (aj — Oj)/y/2i (with j = 1, 2) are the quadratures for 
the two lasing modes. 

After substituting the definitions of u and v in Eq. (SI6), 
we obtain 


It is clear that the equations of motion for the expectation val- (AM) + (AM) 
ues (a\) and {a\) can be obtained by taking the conjugates of 


2[1 + + (a 2 aj) + (aiaj) + (cl^cl^) 

~(ai){af) - ( a 2 )(at ) - (ai)(a 2 ) - (af)<^) 



(Am ) 2 +(Am) 
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We can calculate (A u ) 2 + (Av ) 2 versus time using 
Eqs. (S11)-(S 16), their conjugate equations and the experi¬ 
mental parameters. The data are presented in Fig. SI. We 
thus show the possible existence of quantum correlation be¬ 
tween the two lasing fields in our CEL system. 


We now define the polar coordinates n and 6 i (I = 1,2) via 
relations as in Ref. [1] 

q { = Vi exp (iOi) 

We define the difference and mean angle variables as 
0 = 

0 \ H - 0 2 

V = -O-# 


Phase diffusions 


and then we have 

d exp (—iOi) 
dqi 2 


d Id 
dri + in dOi 


To study the phase noise, we can convert Eq. (S2) into 
an equivalent Fokker-Planck equation for the P-representation 
P = P(qi, q*, q 2l q 2 ) using the following relations: 
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Under the lasing conditions, the variation in the photon num¬ 
ber is small and we can consider that Jr- = 0. Then, we have 

dri ’ 
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characterize the phase diffusions of 0 and 


Here, Dqq and D 
r /, respectively. 

From the Fokker-Planck equation, we can derive the equa¬ 
tions of motion for the phase difference 0 and phase sum r] as 
follows: 


dB_ 

dt 

dr] 

dt 


As shown in Fig. 1(A), (B) and (C) of the main text, the 
single three-level artificial atom with cyclic transitions [3] is 
based on a “tunable gap flux qubifi’circuit [5, 6] capacitively 
coupled to a multimode transmission line resonator (TLR) via 
C c [7]. The effective Josephson and charging energies of each 
junction shown in blue in Fig. 1C are Ej = h x 88 GHz and 
E c = h x 3.1 GHz (where E c = e 2 /2 C for the junction ca¬ 
pacitance C), a = 1.02 (ratio of the total junction capacitance 
marked in red in Fig. 1C to C), C c = 3.1 fF, and the area ratio 
of the main loop to the a-loop is 6.1. The niobium (Nb) TLR 
is patterned by the dry etching of a 50-nm-thick Nb thin film 
deposited on a thermally oxidized silicon substrate. The Nb 
TLR structure has a center conductor width of 10 fi m and gaps 
of 6 fi m, resulting in a wave impedance of about 50 U with a 
exp [i2r)\ + H c t3ta ^ l en glh °f the central conductor of L ~ 10.2 mm. Our 
Experiment is carried out in a dilution refrigerator at a base 
temperature of about 30 mK. 

The artificial atom is fabricated near the voltage antin¬ 
ode of the resonator by the electron-beam lithography and 
exp [—i2r]\ +AI/A Ox/Al shadow evaporation techniques. The energies of 
the three lowest levels of the atom (ground | g), first excited 
state \e) and second excited state | d)) depend on the external 
magnetic flux, reaching minima for the transitions from | g) 
to \e) and | d) at half-integer flux quanta 4>jv = {N + 
(where N is an integer). Compared with superconducting sys¬ 
tems based on the conventional flux qubit geometry [8, 9], our 
atom has additional tunability due to the implementation of an 
a-loop - a dc SQUID, which allows the qubit energy gap A to 
be tuned by controlling the magnetic flux through the SQUID 
iqlgop. When the biased flux is close to the two low- 
estJenergy eigenstates (| g) and |e)) of the artificial atom are 
he superposition of the clockwise and anticlockwise circu- 
ng current states of the main loop. The atomic transition 
frequency in the vicinity of at in the first order is described 
by hu eg 


exp [i2rj\ + H.c. 


H.c 


H.c 


exp (i 2 rj) 


exp [i 2 r]\ + H 


= D e , 

(S28) 

= A,- 

(S29) 


y/( 2 I P S4 >) 2 + A(4>tv) 2 , where 54> = — 4>at and 

I p is the persistent current in the main loop. Here, we neglect 
the weak dependence of A on <f> when 54> <C $0 by assum¬ 
ing that A(4>) « A(<f> 7 v)- By choosing and £4>, we can 
adjust both frequencies, uj eg and cc^ e , and bring them close to 
the resonance of the two lowest modes of the resonator, which 
are two independent fixed-frequency oscillators with the fre¬ 
quencies uj\ and l>j 2 . 


Extracting the dissipation terms of the artificial atom 


Here, Dq = 0 or D g = 0 means that the phase is locked. 
However, Dqq < 0 or < 0 means that the phase diffu¬ 
sions are suppressed. We substitute all the dissipation terms, 

r% = yj(a~i oi) = Vb and r 2 = = \[2 to calculate 

Dpp. Finally, we find that D vri < 0, which means that the 
mutual phase-diffusion noise is suppressed. 


As shown in Fig. S2, we simulate the transmission spectrum 
for the second mode of the resonator around 12 GHz in Fig. 
2B to extract the dissipation terms of the three-level artificial 
atom at the CEL working point. For a probe field with the 
frequency ujd driving the resonator H^ = ^^(aj e~ luJdt — 
a 2 e luJdt ), the total Hamiltonian in the rotating frame at the 
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frequency is 


12.07 


H = h(U2-U d )a%a2 + -{u>dg-CU d )(Jdd+ 7;( UJ eg-U d )(T ee + ',^ 

Z Z N 

a de + &ed) + %2 (^2 G eg + Gge'j + ^5(^2"C) 

> 

where the amplitude is 0,^ = VN in terms of the photon c 

number TV created by the driving field. To extract the dissi- =3 

pation terms of the artificial atom, we can simulate the trans- cu 

mission spectrum in Fig. 2B by solving the stationary master u- 
equation p = 0 with the Hamiltonian in Eq. (S30) and the de¬ 
cay rate of the second mode. The transmission amplitude t in 
the second-harmonic mode of the resonator can be written as 


12.04 

12.01 

11.98 

11.95 


11.92 


t = — 


in 2 (^ 2 ) 

^d 


(S31) 


where ( 02 ) is the expectation value of the photon field in the 
second-harmonic mode of the resonator. In the weak driving 
limit, Eq. (S31) returns back to Eq. (9) in Ref. [4]. We obtain 
the parameters 721 ~ 2tt x 1.5 MHz, 722 = 27r x 6 MHz, 
732 = 27rx3MHz, 731 = 27T x 8 MHz and 733 = 27tx4MHz. 
We consider that the small deviation between the experiment 
and the theoretical simulation is due to the fact that Eq. (S30) 
is not the full system Hamiltonian. 
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Using the sample circuit parameters, we can calculate the 
coupling strength between the different transitions of the atom 
and the different modes of the resonator as shown in Fig. S3. 
In the figure, gi, g 2 , gs, g 4 and g$ correspond to the inter¬ 
actions | 0 d 0 ) <=> |le 0 ), | 0 e 0 ) <=> \0gl), | 0 e 0 ) <=> 

|IgO), |0d0) <=> |0el) and the coupling strength between 
| 0 d 0 ) <=> | 0 g 0 ) and the third mode of the resonator, re¬ 
spectively. At the working point, g\ = 2tt x 90 MHz, g<± = 
27r x 78 MHz and g$ = 2 tt x 225 MHz. 

Experimental setups for the transmission and emission 
measurements 


FIG. S2: Experimentally measured (A) and theoretical simulations 
(B) of the transmission spectrum for the second mode. From the 
simulations, we can extract the dissipation terms of the three-level 
artificial atom. 


in Fig. S7 A. 

The experimental setup used to demonstrate the correlation 
of the lasing fields with two diffent colors in the two-mode 
CEL is presented in Fig. S7 B. It is necessary to seperately 
balance the gains in the two different chains to obtain maximal 
modulation in the interference fringes. 


The experimental setups used to measure the transmission 
and emission are presented in Fig. S4 A and Fig. S4B, respec¬ 
tively. 

In the experiment, when we change the pumping power 
and pumping frequency, the sum of the two center frequen¬ 
cies of the emission spectra is always equal to the pumping 
frequency, as shown in Fig. S5 and Fig. S 6 . 

Experimental setup for measuring the quenching of the phase 
diffusion and the correlation of the two-mode CEL 

The experimental setup used to demonstrate the quenching 
of the phase-diffusion noise of the two-mode CEL is presented 
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FIG. S3: Theoretical calculation of the coupling strength between 
different atomic transitions and cavity modes, gi, < 72 , g?>, g 4 and 
#5 correspond to the interactions |0d0) <=> |leO), |0e0) <=> 
\0gl), |0e0) <=> \lgO), |0d0) <=> |0el) and the coupling 
strength between |0d0) <=> |0g0) and the third mode of the res¬ 
onator, respectively. 
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FIG. S4: (A) The experimental setup for measuring transmission 
spectrum and energy spectroscopy of the coupled system. (B) The 
experimental setup for measuring two-mode CEL emission spec¬ 
trum. 
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FIG. S5: Sum of the two center frequencies of emission spectra is 
equal to the pumping frequency with changing the pumping power. 
Emission spectra from the first mode (A) and second mode (B) with 
changing the pumping power. (C) Sum of the two center frequencies 
of the emission spectra in (A) and (B) with changing the pumping 
power. 
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FIG. S6: Sum of the two center frequencies of emission spectra 
is equal to the pumping frequency with changing the pumping fre¬ 
quency. The emission spectra from the first mode (A) and the second 
mode (B) with changing the pumping frequency. (C) The sum of 
the two center frequencies of emission spectra in (A) and (B) with 
changing the pumping frequency. 




FIG. S7: (A) Experimental setup for measuring the quenching of the 
phase-diffusion noise of the two-mode CEL. (B) Experimental setup 
for measuring correlation of the two-mode CEL. 
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